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LARGE  VISCOPLASTIC  DEFLECTIONS  OF  IMPULSIVELY 
LOADED  PLANE  FRAMES1 

2 3 

by  P.  S.  Symonds  and  C.  T.  Chon 


■anuniN/mtusain  cc?Et 


Abstract 

Applications  are  described  of  two  estimation  techniques  to  obtain  final 
deflections  and  response  times  of  plane  rectangular  frames  subjected  to  im- 
pulsive loading  on  the  transverse  (beam)  member.  Deflections  up  to  roughly 
one  third  the  span  (thirty  thicknesses)  are  estimated  by  the  mode  approxima- 
tion and  deflection  bounds  techniques , treating  the  plastic  rate  dependence 
by  means  of  homogeneous  viscous  constitutive  equations.  Comparisons  are  made 
with  recent  test  results,  and  the  degree  of  agreement  is  discussed  in  terms  of 
the  known  error  sources  of  the  two  techniques. 
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LARGE  VISCOPLASTIC  DEFLECTIONS  OF  IMPULSIVELY  LOADED  PLANE  FRAMES 

1.  Introduction 

The  impulsively  loaded  plane  frame  considered  here  is  assumed  to  exhibit 
strong  plastic  rate  sensitivity  and  to  reach  large  deflections.  We  shall  des- 
cribe applications  of  the  mode  approximation  technique  and  of  the  deflection 
bound  method  to  this  type  of  structure.  Experiments  intended  to  check  on  these 
methods  have  recently  been  made  [1],  The  comparison  between  predictions  of  the 
estimation  techniques  and  the  test  results  is  discussed  in  order  to  assess  the 
relative  importance  of  the  intrinsic  errors  of  the  methods,  and  those  due  to 
further  idealizations  and  approximations  made  in  their  application  to  this  type 
of  structure. 

The  two  types  of  frames  considered  are  shown  in  Fig.  1.  Type  (a)  in  Fig. 
la  has  a concentrated  impulse  applied  to  a small  block  at  the  midpoint  of  the 
beam  (transverse)  member,  while  type  (b)  in  Fig.  lb  has  a distributed  impulse 
applied  over  this  member.  These  loads  are  idealized  as  impulsive  (zero  duration), 
imparting  specified  initial  velocities  with  negligible  initial  displacements. 
Symmetric  deformations  are  expected,  the  main  displacement  magnitude  being  the 
displacement  at  the  midpoint  C of  the  beam  member.  At  large  deflections  such 
as  indicated  in  Fig.  2b,  a finite  lateral  displacement  occurs  at  the  top  B of 
each  column.  The  particular  frames  studied  here  and  in  the  related  tests  had 
lengths  L^  = 5.625  in.  in  the  case  of  the  "type  (a)"  frames,  L^  = 6.00  in.  for 
the  type  (b)  frames,  with  L^  = 8.00  in.  in  both  cases.  Two  strongly  rate- 
sensitive  metals,  mild  steel  and  commercially  pure  titanium  were  used.  In  these 
circumstances  the  large  deflections  have  primarily  kinematic  effects,  without 
requiring  drastic  changes  in  the  stress  field;  the  response  remains  primarily 
flexural. 


2.  Basic  Concepts  and  Equations 

The  essential  concepts  of  the  two  techniques  will  be  outlined  briefly.  Fur- 
ther information  on  the  approaches  may  be  found  for  the  bound  method  in  [2,  3] 
and  for  the  extended  mode  apnroximation  technique  in  [3,  4], 

The  deflection  bound  method  for  a structure  of  viscoDlastic  material,  sub- 
jected to  impulsive  pressures  at  t $ 0 and  thereafter  unloaded t requires  the 
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solution  for  quasi-static  deflections  due  to  a certain  concentrated  force  P^n  : 
if  the  deflection  is  sought  at  point  A of  the  structure  (located  by  appropr- 

«s 

iate  coordinates  in  direction  n (a  unit  vector),  is  applied  at 

« sf 

this  point  and  in  this  direction.  If  the  force  P^n  is  such  that  it  does  work 
at  least  ecual  to  the  specified  initial  kinetic  energy  of  the  structure,  then 
the  displacement  on  which  it  does  work  is  an  upper  bound  on  the  same  deflection 
quantity  in  the  response  to  the  given  dynamic  loading.  If  t^.  is  the  time  at 
which  the  plastic  deformation  is  completed,  the  theorem  requires  the  work  of  the 

g 

force  P^  to  be  done  in  time  tf  . ( For  a time-dependent  material  the  work  to 
reach  specified  final  strains  depends  on  the  time  in  which  they  are  reached.) 

In  mathematical  terms  the  bound  theorem  states  that 
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where  u^  is  the  final  displacement  at  point  A in  direction  n of  the  dy- 
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namically  loaded  structure,  reached  at  time  tf  ; and  uftn  the  displacement  at 


time  t,  due  to  the  static  force  P7 
f An 


the  inequality  holds  provided 
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where  Kq  is  the  initial  kinetic  energy  of  the  impulsively  loaded  structure 

S 

and  the  right-hand  side  represents  the  work  done  by  the  force  P^n  in  the  in- 

y g f 

terval  of  time  0 < t < t^  ; W(q^  ) denotes  work  per  unit  volume  written  as 
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function  of  the  terminal  generalized  strains  q^  » 1 = l,..r,  for  general- 
ized strain  and  stress  states  having  r components. 

As  noted  above,  this  work  depends  on  the  interval  t^  . Generally  it 
also  depends  on  the  path,  i.e.,  on  the  sequence  of  strain  states  from  the  in- 
itial to  the  terminal.  We  have  eliminated  path  dependence  by  using  concepts 

V 

of  minimum  work  naths  [2,  3],  W denoting  the  work  per  unit  volume  which  is 
a minimum  for  given  terminal  strain.  The  total  work  is  evaluated  as 
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where  T(Q^  ) is  a homogeneous  function  of  terminal  stresses  0^  , whose  de- 

gree of  homogeneity  is  n + 1.  More  generally,  the  function  V furnishes  con- 
stitutive equations  relating  stress  to  strain  rate  states  in  a convenient  form, 
from  the  property 


at 

q1  " 30, 


(2) 


where  q^  are  strain  rate  components. 

The  "intrinsic  error"  of  the  deflection  bound  is  positive:  an  upper  bound 
is  obtained.  In  this  statement  the  comparison  is  made  between  a computed  quasi- 
static deflection  according  to  a certain  mathematical  model,  and  a deflection 
resulting  from  impulsive  loading  on  the  structure  represented  by  the  same  model. 
When  comparison  is  made  instead  with  the  final  deflection  of  an  actual  structure 

l 

as  opposed  to  a mathematical  model,  the  idealizations  and  approximations  adopted 

sf 

in  the  calculation  of  uAn  may  introduce  further  errors.  Some  of  the  ideali- 
zations are  indicated  above;  they  will  be  discussed  with  others  more  fully  in 
a later  section.  Obviously  the  aim  is  to  make  the  calculation  in  such  a way 
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that  the  incidental  errors  are  minor  compared  to  the  intrinsic  ones.  The  need 
for  comparisons  with  experiments  rather  than  merely  with  computer  outputs  is 
also  obvious. 

The  essential  concept  of  the  mode  technique  is  that  of  obtaining  an  ap- 
proximation to  the  actual  response  from  a simpler  solution  which  satisfies  all 
the  field  equations  (dynamics,  kinematics  including  boundary  fixing  conditions, 
and  constitutive  equations),  but  disagrees  in  general  with  the  imposed  initial 
velocities.  Such  simpler  solutions  can  be  found  under  certain  conditions  in 
modal  form,  with  velocity  field,  for  example,  written  as 


u?(x,t)  = T(t)$.(x) 


(3) 


where  i * 1,  2,  3;  T(t)  is  a scalar  function  of  time  t;  and  ♦j(x)  is  a 
vector-valued  function  of  space  coordinates  x . The  initial  velocities 


u^(x,0)  = T^tx 


) 


(4) 


where  Tq  » T(0),  usually  differ  from  the  velocities  specified  at  t = 0,  namely 


u-^x.O)  = u?(x) 


(5) 


since  the  shape  functions  ♦^(x)  are  properties  of  the  structure.  However 


the  initial  magnitude  T can  be  chosen  optimally  by  taking  [5] 

o 
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(6) 


where  p is  the  mass  density  and  the  integral  is  over  the  volume  V of  the 

structure.  This  value  of  T minimizes  the  initial  magnitude  of  the  follow- 

o 

ing  functional: 

A(t ) = A[Uj(x,t ) - u.(x,t)]  * p(Uj  - UjMu.  - Uj)  dV  (7) 
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It  has  been  shown  [6]  for  a wide  class  of  material  behavior  of  essentially 
viscous  type  (generalized  strain  rates  written  as  functions  of  generalized 
stress)  that  A(t)  is  a non-increasing  function,  and  decreases  whenever  plastic 
flow  occurs  and  the  two  solutions  do  not  have  identical  stresses  or  strain  rates. 
Thus  the  actual  velocity  field  and  that  of  the  mode  solution  approach  each  other 
in  this  sense.  With  Tq  chosen  according  to  Eq.  (6),  the  two  solutions  may 
become  identical  after  a certain  interval.  The  final  major  displacement  com- 
puted from  the  mode  solution  is  therefore  usually  much  closer  to  that  of  the 
structure  than  the  initial  mode  amplitude  io  to  the  corresponding  given  initial 
velocity. 

The  fined  major  deflection  of  the  mode  response  will  be  greater  or  less 
than  that  of  the  structure  whose  initial  velocity  is  the  specified  one,  depend- 
ing on  whether  the  initial  mode  velocity  field  is  more  or  less  concentrated  than 
the  given  velocity  distribution.  This  intrinsic  error  evidently  can  be  posi- 
tive or  negative.  Here  again  the  "error"  refers  to  quantities  computed  using 
a certain  mathematical  model  to  represent  a structure.  The  calculation  of  the 
mode  response  involves  idealizations  and  approximations,  and  when  the  results 
are  comp  aired  with  deflections  and  response  times  of  a real  structure*  further 
errors  may  appear.  These  will  be  discussed  in  a later  section. 

The  mode  method  requires  the  integration  of  the  equations  of  the  structure, 
from  the  initial  mode  form  velocity  field  to  the  end  of  the  motion.  This  in- 
tegration presents  no  difficulty  if  the  equations  allow  the  separated-variable 
form  of  Eq.  (3)  to  hold  during  the  entire  response.  If  small-deflection  eq- 
uations are  not  used,  such  "permanent"  mode  solutions  do  not  exist.  A conven- 
ient way  of  performing  the  integration  in  these  cases  is  by  means  of  a sequence 
of  "instantaneous"  mode  form  solutions  [3,  4].  At  finite  deflections  the  so- 
lution can  be  put  in  the  form  of  Eq.  (3)  if  the  deflection  is  regarded  as  in- 
stantaneously fixed  and  known.  Thus  the  instantaneous  shape  function  $.(x) 
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<t> 

changes  during  the  response  and  may  be  written  as  ^(x).  Successive  "solutions” 
in  this  fjro  are  linked  through  equations  of  the  form 

a*  <tn>  <tn+l> 

ui*x,tn+l*  3 ui*x,tn*  + T~  C^t^^Cx)  + T^tn+l^*i(x^  (8a) 

T<Vi>  ' T<V  * T-  [*<V  * *<W]  <8b> 

This  method  is  not  exact,  even  though  the  field  equations  are  instantaneously 
satisfied;  the  mode  solution  has  a smaller  energy  dissipation  rate  than  the  actual 
motion  at  the  same  level  of  kinetic  energy  [7].  This  method  may  be  expected  to 
lead  to  a response  of  longer  duration  and  larger  deflections  than  the  actual 
motion  of  the  structure  would  be  under  the  same  starting  conditions,  so  that  the 
approximation  due  to  this  device  for  carrying  out  the  integration  for  large  de- 
flections is  therefore  expected  to  be  conservative. 

Here  we  write  the  equations  used  in  both  the  mode  technique  and  the  deflec- 
tion bound  method;  the  dynamical  equations  for  the  latter  are  obtained  by  setting 
the  accelerations  equal  to  zero. 

The  notation  is  indicated  in  Fig.  2b.  We  use  rectangular  coordinates  x,  y 
with  origin  at  the  base  A of  the  left-hand  coluim.  Velocity  components  are 
shown  in  Fig.  2b.  Nondimensional  quantities  are  used,  defined  as  follows 


x 

X*L, 


w 

w 3 H 


u = L aji 


W - — — 


where  x,  y;  t;  u,  w are  coordinates,  time,  and  displacements  in  physical 
units,  H is  thickness,  and  t = p/oQ  is  a reference  time,  p being 
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mass  densitv  and  oq  beinp  a stress  property  obtained  from  tests  on  strain 
rate  dependent  plastic  behavior,  as  defined  more  fully  below.  For  general 
stress  and  strain  states  in  a one-dimensional  structure,  relations  between 
bending  moment,  axial  force  and  the  corresponding  strain  rates  are  required. 
We  take  the  viscoplastic  behavior  to  be  expressible  in  the  form 


e > 0 


(10) 


where  o , e are  uniaxial  stress  and  plastic  strain  rate,  respectively,  cor- 
responding (in  general)  to  a fixed  level  of  plastic  strain  ; and  oq  , 
eQ  , n are  experimental  constants  appropriate  to  that  level.  For  mild  steel 
the  pair  o , e are  more  appropriately  taken  as  lower  yield  stress  and  cor- 
responding strain  rate,  respectively;  then  oq  has  the  significance  of  yield 
stress  at  zero  3train  rate.  Equation  (10)  is  capable  of  very  good  represent- 
ation of  observed  dynamic  plastic  behavior  for  the  metals  with  strong  strain 
rate  sensitivity,  provided  strain  rate  history  effects  are  negligible.  General- 
izations of  Eq.  (10)  are  easily  written,  e.g.  [4,  8].  These  involve  a yield 
condition,  plastic  strain  rates  being  zero  for  stress  states  inside  a yield 
surface.  We  adopt  a simpler  form  derived  from  that  of  Eq.  (10),  which  is  homo- 
geneous and  involves  no  yield  condition  [9].  This  "matched  viscous"  represent- 
ation replaces  Eq.  (10)  by 


where  o'  = pa  , n*  = vn  , with  (lib) 

o o 


e 

v = — 
o 

e 


(lid) 


where  e is  the  strain  rate  at  which  the  two  expressions,  Eq,  (10)  and  (11a), 
are  matched  in  the  sense  of  having  the  same  stress  and  slope  do/de.  In  gener- 
alized  forms,  c is  the  appropriate  effective  strain  rate. 

The  following  expressions  are  generalizations  of  Eq.  (11a)  appropriate  for 
general  states  of  stress  in  our  frame: 


1 f.-  -,1/n'  ,*  >• 

m = j-vm  + n|  sgn(C  + n)  + K 

s = + sgn( £ + n)  - |C 


. 1/n'  . . ) 
n|  sgn (l  - n) f 

1/n'  . . 

n|  sgn(C  - n)> 


where  m = 


S = N • 


M a —■  a , N = bHo  , 
o 4 o o o 


k = u e 
o H o 


(12a) 


(12b) 


(12c) 


(12d) 


M , N being  physical  bending  moment  and  axial  force,  respectively  and  ic  , e 
being  curvature  rate  and  axial  strain  rate  respectively  in  physical  terms.  Al- 
though derived  from  a sandwich  beam  model  [2,  3,  4],  the  values  for  pure  bend- 
ing and  extension  are  either  exact  or  extremely  close  to  the  correct  values  for 
a rectangular  solid  section.  These  constitutive  equations  are  conservative  in 
the  sense  that  for  a given  strain  rate  state  the  stress  levels  are  less  than 
those  in  the  solid  section. 

We  next  write  the  strain  rate  components  in  terms  of  velocity  and  displace- 
ment components  in  the  two  members  AB  and  BC,  using  dimensionless  variables 
(see  Eqs.  (9,  12)  and  Fig.  2): 
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5AB  " a 2 UAB 
L2 


• 1 • 

r s - — m" 

^BC  o BC 


A -I  if  Si'  ♦ u'  u' 

nAB  " a 72  [“  AB  ABuAB 


(13a, b) 


(13c) 


V = ? (7i  “bc  * wbc"bc] 


where  o = 4L2e  t/H2  = (8L^e  /H2)/o/o 
1 o 1 o o 


(13d) 


and  a prime  denotes  differentiation  with  respect  to  y in  AB  and  to  x in 
BC;  note  that  0 f y f 1 , 0 < x < 1 . The  expressions  for  axial  strain  rate 

in  Eqs.  (13c,  13d)  contain  terms  where  the  deflection  curve  appears  through  the 
rotation  (slope),  as  in  the  "von  Karman  equations"  of  buckling  theory.  They 

are  second  order  correction  terms,  valid  for  deflections  of  moderately  small  size 
compared  to  the  span. 

The  equation  of  energy  -dissipation  rate  can  be  written  for  both  types  of 


frame  as 


l r1  r1 

- ET  (“ab5ab  + WAB)dy  - (VV  + VV)dx  - kVl 
-T)  'o 


= ° L~  (mAB5AB  + SABnAB)dy  +“  (mBC5BC+  SBCnBC)dX  (14) 

1 -'0  '0 

where  k = G/2pbHL^,  and  w.^  = wfiC(l,t),  w^=  w"BC(l,t)  are  velocity  and  accel- 
eration of  the  block  of  mass  G.  When  applied  to  the  frame  with  attached  mass, 
L is  given  the  reduced  value  5.625  in.;  for  the  frame  without  attached  mass, 
k = 0 and  * 6.00  in. 

The  equations  of  dynamics,  end  conditions,  and  equations  of  continuity 
consistent  both  with  the  kinematic  forms  of  Eqs.  (13)  and  the  dissipation- 


energy  rate  equation  (14)  are  derived  as  the  Euler  equations  of  the  latter, 

the  velocities  and  associated  strain  rates  being  treated  as  virtual  quantities. 

The  dynamical  equations  are  found  to  be 

,2 


mAb  + 4(sABUAB) ' = 72  UAB 
L1 

mK  + 4(sbcwbcv  = V 

L2  , L2  .. 

H AB  2 AB 

b , . 1 

4 H S2C  " UBC 


(15a) 

(15b) 

(15c) 

(15d) 
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End  conditions  are: 


At  C,  x = 1 : 


mBC ( 1 *t ) = “ k”BC(1,t) 


(16a) 


BC 


(l,t)  = w^c(l,t)  = Upp( 1 ,t ) = 0 


BC 


(16b,  c,  d) 


Equation  (16)  is  recognized  as  the  equation  of  motion  of  the  attached  block; 
k is  put  equal  to  zero  for  the  frame  without  attached  mass. 


At  A,  y = 0 : uAR(0,t)  = u^0(o,t)  = wftB(0,t)  = 0 


At  B,  y = 1,  x = 0;  continuity  conditions  are  found  to  be 


UA3(1’t)  = UBC(°’t);  WAB(1’°  = WBC(°’t) 


mAB^1,t^  = mBC^0,t^’  L2  UAB^1,t:^  = WBC^0,t^ 

I. 

* "•Ab(1  '’"iB0’'1  ' V SBC(0t)  “ ° 


(17a,  b,  c) 


(18a, b) 
(18c, d) 
(18e) 


.... 
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We  have  made  the  assumption,  in  calculations  made  to  date,  that  the  deform- 


ations are  flexural,  with  unimportant  effects  due  to  lengthening  or  shortening 


of  frame  members.  This  is  reasonable  in  view  of  the  proportions.  Hence  the 


simplifying  assumption  of  inextensibility  was  adopted,  with  n = 0 in  both 


members.  The  axial  forces  are  then  reactions,  and  Eqs.  (12c,  d)  furnish  re 


The  constitutive  equations  then  reduce  to  those  for  pure  bending 


In  the  mode  approximation  technique  we  look  for  solutions  in  separated 


variable  form,  writing 


w „(l,t)  is  the  transverse  velocity  at 


1 at  x = 1:  thus  w. 


midpoint  C (the  major  velocity  magnitude).  With  these  forms  substituted  in  the 
equations  of  dynamics  Eqs.  (15)  and  the  kinematic  relations  Eqs.  (13),  it  may  be 


seen  that  the  system  of  equations  can  be  split  into  ordinary  differential  equa 


tions  containing  either  space  or  time  variables,  provided  the  rotation  terms  from 


known  at  some  stage  of  the  response,  the  resulting  eigen-problem  can  be  solved 


to  furnish  the  corresponding  shape  functions  $ and 


and  the  current 


acceleration  w.,{  and  velocity  w, 
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Details  of  the  integration  by  an  iterative  scheme  are  given  in  the  Appendix. 

The  deflection  bound  is  determined  by  a closely  related  numerical  scheme,  also 
outlined  in  the  Appendix. 

The  integration  of  the  nondimens ional  equations,  starting  from  an  initial 
velocity  field  w^$(x),  can  be  completed  once  values  are  assigned  of  the  par- 
ameters a , n , L^/L ^ , L^/H  , as  well  as  of  the  initial  velocity  amplitude 
w^  . (Note  that  u and  n'  = vn  are  obtained  from  these  and  the  current 
strain  rates).  The  end  of  the  motion  occurs  at  time  t*  such  that  w^(t*)  = 0 , 
and  the  final  deflection  components  are  w^  , w^2  . The  final  deflection 
amplitude  and  response  time  can  be  presented  as 

-f 

wft  f 

— wA  = * 0 » n , L^/H , L^/Lj  » K)  (21a) 

t* 

— =t*=F2(w*,  a,  n,  Lj/H,  , k)  (21b) 

where  we  recall  t = 21^/  p/oq  , a = (8L^eq/H2) J p /oq  , fc  = G/2L^pbH. 

Examples  of  plots  of  w^  and  t*  as  functions  of  w£  for  various  values 

of  the  material  and  geometrical  parameters  are  shown  in  Figs.  3 and  4.  The  de- 
crease of  velocity  amplitude  to  zero  is  illustrated  in  Fig.  5.  The  two  sets 
of  curves  of  Fig.  3 illustrate  how  the  parameters  a and  n affect  the 
final  deflection,  the  remaining  parameters  of  Eq.  21  being  held  fixed.  The 
curves  of  Fig.  4 indicate  how  the  duration  time  is  affected.  Weak  dependence 
on  both  a and  n is  evident.  For  example,  multiplying  o by  1/4  or  1/100 
reduces  the  deflection  by  only  about  10  percent  or  50  percent,  respectively. 

Fig.  3b  shows  that  dependence  on  n is  similarly  weak.  The  insensitivity 
to  these  parameters  can  be  understood  from  the  form  of  the  equations.  For 
example,  the  energy  rate  - dissipation  equation  based  on  the  mode-fom 
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where  the  line  integrals  extend  over  the  half-frame  ABC  , and 
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Here,  as  in  all  equations  of  the  system,  a appears  only  in  the  combination 
(wA/a)^n  . Since  n'  = vn  , where  v depends  on  current  strain  rate  (Eqs. 
(11c))  and  is  often  approximately  2;  and  since  n is  large  (5  for  steel  or  9 for 
titanium),  the  weak  dependence  on  a is  evident.  We  obtain  results  for  a 
perfectly  plastic  material  by  putting  y = 1 , n'  ; as  expected,  then  the 
acceleration  is  independent  of  velocity  and  the  parameter  a disappears. 

This  insensitivity  to  a and  n is  an  advantage  in  applications.  Once 
the  integration  has  been  carried  out  for  nominal  values  of  the  parameters , giving 
curves  of  w£  and  t^  as  function  of  w^  over  a suitable  range,  these  curves 
can  be  used  for  many  particular  cases  of  interest.  For  the  uniform  frame  two  such 
curves  are  shown  in  Fig.  6 for  "steel"  and  "titanium",  and  an  analogous  pair  of 
curves  is  shown  in  Fig.  7 for  the  frame  with  attached  mass.  An  inverse  procedure 
is  used  to  apply  these  to  the  data  of  a particular  test.  First,  the  initial 
velocity  matching  formula  Eq.  (6)  of  the  mode  approximation  is  written(in  terms 
of  nondimens ional  quantities).  Then  from  the  relation  between  measured  impulse 
I in  physical  units  and  the  non-dimensional  initial  velocity  wc  , the  impulse  I 
corresponding  to  a particular  w?  is  obtained.  Thus  we  obtain  w£/H  and 
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t * tfA  as  functions  of  I.  The  results  for  the  two  types  of  frames  are 
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Taking  data  for  the  two  types  of  frame  and  two  materials,  final  deflection 
ratios  are  shown  in  Figs.  6 and  7 as  functions  of  nondimensional  mode  velocity 
w^  . From  these  using  Eqs.  (23b,  d)  are  obtained  curves  of  deflection  versus 
impulse  in  physical  units.  Examples  of  such  curves  are  given  in  Figs.  8-12,  with 
points  shown  also  from  tests  [1]. 

In  computing  the  Impulse  I corresponding  to  a chosen  value  of  wj  from 

Eq.  (23b ,d),  values  of  b,  H,  p and  oq  appropriate  for  a particular  test  series 

are  used,  and  these  are  used  also  to  determine  the  appropriate  value  of  the  par- 
3 • 2 . 

ameter  o = (8L,e  /H  )/  p/o  . However  in  view  of  the  insensitivity  of  the  curves 
10  o 

for  w^  versus  w^  to  a , it  is  evident  that  a curve  for  a nominal  value  of 
a can  be  used  for  other  cases.  Thus  the  curves  of  Figs.  6 and  7 serve  as 
"master  response  curves"  for  frames  of  a certain  type  of  material  and  geometrical 
configuration . 

To  illustrate,  a nominal  value  of  a for  the  steel  frames  tested  [1]  may 
be  taken  as  550,  corresponding  to  material  and  geometrical  values  in  Table  1 of 
[1].  In  particular,  oq  = 33,100  psi  was  used,  from  strain  rate  data  on  the 
material  of  those  tests.  Thus  for  example,  a final  deflection  ratio  wj^  = w^/H  = 20 
corresponds  to  a nondimensional  initial  mode  velocity  w£  = 18.0.  Using  Eqs. 

(23)  with  H = 0.123  in.,  a final  deflection  of  2.46  in.  is  estimated  to  require 
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a uniformly  distributed  impulse  I = 0.75  lb-sec,  as  plotted  in  Fig.  10.  Now 
suppose  a frame  of  the  same  dimensions  but  of  higher  strength  steel  with  oq  = 
90,000  psi  is  to  be  considered.  The  new  a would  be  550/^3.1/90  * 334,  but 
the  same  "master  curve"  of  Fig.  6 can  be  used,  since  the  change  due  to  a 40 
percent  reduction  in  a is  negligible,  by  Fig.  3(a).  Hence  the  impulse  on 
the  stronger  frame  required  to  produce  the  same  deflection  can  be  estimated  as 
0.75/90/33.1  = 1.2  lb-sec. 

This  simple  type  of  calculation  was  used  in  analyzing  the  tests  [13  to 

estimate  the  effect  of  strain  hardening  in  the  titanium  frames.  The  values 

• • 
of  oq  , nf  and  cq  depend  on  the  plastic  strain  at  which  (a  , c)  are  ob- 

served  in  tests  at  nominally  constant  strain  rate.  At  e*3 * 5  = 1 percent  and  2 
percent,  values  of  oq  were  determined  as  35,000  psi  and  38,000  psi,  respect- 
ively. Values  of  total  impulse  for  both  values  are  shown  for  the  estimated 
deflection  curves  for  titanium.  Figs.  9 and  11.  In  most  applications  one  sim- 
ilarly wants  to  explore  a range  of  values  of  various  parameters,  including  var- 
ious load  distributions  over  the  structure.  It  is  seen  that  the  present  tech- 
nique is  particularly  efficient  for  such  explorations,  and  contrasts  with  the 
typical  wholly  numerical  approach. 

3.  Comparisons  with  Test  Results 

As  already  noted,  experimental  results  [1]  are  now  available  for  check- 

ing the  estimation  techniques.  The  basic  concepts  (minimum  potential  energy 
and  the  mode  velocity  matching  device  of  Eq.  (6))  are  not  in  doubt,  but  the 
magnitudes  of  the  intrinsic  errors  in  estimates  of  deflections  and  response 
times  are  unknown.  Unless  these  errors  are  reasonably  small,  and  unless  the 
further  errors  involved  in  implementing  the  techniques  are  also  reasonably 
small  as  well  as  generally  predictable  as  to  sign,  the  methods  will  remain 
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of  uncertain  meaning  and  limited  usefulness.  Thus  the  experiments  are  crucial 
for  assessing  the  relative  importance  of  the  two  classes  of  errors , and  for 
confirming  or  otherwise  our  expectations  as  to  the  incidental  sources  of  er- 

* - t \ DJTJ  || 

ror. 

The  deflection  upper  bound  curves  plotted  in  Figs.  8-11  in  every  case  lie 

; ' 1 

above  the  test  points,  indicating  that  the  incidental  errors  in  this  calcul- 
ation are  either  small  or  predominately  positive.  However,  for  the  mode  tech- 
nique the  situation  is  not  so  clear.  For  the  concentrated  impulse  tests  (frame 
type  (a))  the  intrinsic  error  is  negative,  since  w£  < w°  from  Eq.  (23a);  for 

the  uniform  impulse  tests  (on  frames  of  type  (b))  w*  > w°  from  Eq.  (23c), 

c 

. , . • •-  * ' • ■ ■' 

and  the  intrinsic  error  is  positive.  If  all  other  errors  were  negligible,  we 

would  expect  the  tests  points  to  lie  above  the  estimated  curves  (for  finite 
deflections)  in  Figs.  8 and  9,  and  to  fall  below  them  in  Figs.  10  and  11.  No 
such  consistent  relations  are  discernible  in  these  figures.  In  the  case  of  the  steel 
frames  with  either  concentrated  or  distributed  loading  the  test  points  lie  es- 
sentially  on  the  estimated  curve  (although  for  concentrated  impulses  they  tend 
to  be  consistently  higher  at  the  larger  impulse  magnitudes).  For  the  titanium 
frames  under  concentrated  impulses  the  test  deflections  fall  well  below  the  es- 
timated curves,  although  approaching  them  at  the  larger  impulse  magnitudes. 

Under  distributed  impulses  the  titanium  frames  showed  final  deflections  agree- 
ing  quite  closely  with  the  estimated  curves.  (For  titanium  two  estimated  de- 
flection curves  by  the  mode  method  are  drawn,  for  two  values  of  oq  ; as  already 
noted,  these  correspond  to  two  choices  of  strain  level) . The  fact  that  the 
experiments  do  not  clearly  show  the  intrinsic  errors  suggests  that  they  are 
masked  by  effects  of  the  further  idealizations  and  approximations  made  in  apply- 
ing the  mode  technique.  We  next  list  these  and  consider  them  briefly. 
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(1)  A "rigid-viscoplastic"  theory  was  used:  elastic  deformations  were  neg- 
lected (elastic  moduli  taken  as  infinite).  Here  we  are  considering 
strongly  rate  sensitive  metals  so  that  stress  levels  are  raised  and  elas- 
tic strains  probably  increased  in  importance.  On  the  simplest  basis , we 
may  suppose  that  the  energy  ratio  R of  initial  kinetic  energy  to  total 
elastic  strain  energy  capacity  must  exceed  about  6 for  an  error  of  less 
than  15  percent  in  the  permanent  deflection;  this  takes  over  a result  [10] 
derived  from  a one-degree-of-freedom  model  and  seems  conservative  for  per- 
fectly plastic  behavior.  The  values  of  impulse  required  for  R ) 6 are 
indicated  in  Figs.  8-11.  If  the  same  results  held  for  the  present  rate 
sensitive  materials  there  would  appear  to  be  some  correlation  between  this 
(arbitrary)  criterion  for  the  neglect  of  elastic  deformations  an^  the  ex- 
perimental discrepancies.  In  the  one  case  where  the  criterion  is  not 
satisfied  in  some  of  the  tests,  namely  the  concentrated  impulse  tests  on 
titanium  frames  shown  in  Fig.  14,  the  test  points  at  the  smaller  impulse 
magnitudes  fall  considerably  below  the  estimated  deflection  curves,  but 
approach  them  as  the  impulse  is  increased.  Mote  that  for  this  case  the 
intrinsic  error  is  negative  (the  test  points  should  lie  above  the  esti- 
mated curve  in  the  absence  of  other  errors);  hence  the  effect  of  elastic 
deflections  would  be  even  larger  than  indicated. 

As  already  noted,  errors  due  to  elastic  effects  may  be  larger  when  the 
material  is  strongly  rate  sensitive  than  indicated  by  the  criterion  ap- 
plied, as  above,  for  perfectly  plastic  behavior.  The  omission  of  elastic 
deformations  makes  for  substantial  simplification,  but  further  basic  re- 
search is  needed  to  determine  what  errors  are  likely  to  be  caused  by  it. 

It  should  probably  be  regarded  here  as  a major  suspect. 
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(2)  The  constitutive  equations  express  plastic  strain  rates  as  explicit  func- 
tions of  stresses t with  implicit  account  of  strain  hardening  through  the 
experimental  constants.  If  the  arbitrarily  chosen  plastic  strain  level 
is  low  (or  if  the  lower  yield  stress  is  used  rather  than  a fixed  strain 
level),  equations  of  this  type  would  underestimate  strain  hardening  and 
hence  lead  to  over-estimates  of  deflections.  However,  as  illustrated  for 
the  titanium  frames,  it  is  easy  to  deduce  deflection  estimates  based  on 

a larger  magnitude  of  plastic  strain,  and  thus  to  take  account  of  strain 
hardening. 

(3)  Strain  rate  history  effects  were  neglected,  the  experimental  constants 
being  obtained  from  tests  at  nominally  constant  strain  rates.  The  in- 
fluence of  prior  strain  rate  history  has  been  studied  mainly  by  tests  in 
which  the  strain  rate  is  rapidly  increased  [11].  A few  tests  involving 
a rapid  decrease  have  been  described  [12];  these  would  be  more  directly 
relevant  to  the  present  response  histories,  where  the  structure  is  set 
in  motion  in  a very  short  time,  after  which  the  strain  rates  decrease  to 
small  values,  perhaps  roughly  monotonically.  On  the  basis  of  various 
hypotheses  to  account  for  observed  history  effects,  in  particular  the  use 
of  plastic  work  as  a state  variable  [13]  and  dynamic  recovery  [14],  one 
would  expect  strain  rate  history  effects  in  the  present  case  to  lead  to 
positive  errors,  i.e.,  to  cause  the  test  specimens  to  deflect  less  than 
the  estimated  value.  The  data  are  quite  limited,  but  there  seems  little 
reason  to  expect  that  such  errors  would  be  important  here. 

(4)  Constitutive  equations  of  homogeneous  viscous  type  (without  yield  con- 
dition) were  used,  derived  from  inhomogeneous  forms  based  directly  on 
strain  rate  test  data.  The  matching  formula  used  (Eqs.  11c)  is  such  that 
the  stress  levels  are  lower,  at  the  same  strain  rate,  than  those  accord- 


ing  to  the  more  realistic  forms.  Hence  over-estimates  of  deflections  are 
again  expected.  The  errors  due  to  this  device  have  been  very  small  in  the 
examples  worked  to  date  [3*  9].  Experience  is  limited,  but  there  seems 
no  reason  to  expect  large  errors  in  the  present  cases. 

(5)  The  deformation  is  assumed  to  remain  flexural  at  large  deflections,  with 
axial  forces  treated  as  reactions.  Expressions  for  center-line  strain 
rate  in  terms  of  transverse  and  axial  velocity  components  are  of  second 
order  accuracy  (product  of  rotation  by  rotation  rate).  Setting  the  axial 
strain  rate  equal  to  zero  furnishes  relations  between  the  velocity  com- 
ponents, involving  the  transverse  deflection. 

The  errors  due  to  the  use  of  approximate  strain  expressions  at  fin- 
ite displacements,  and  to  the  assumption  of  zero  net  center-line  strain 
rate  are  difficult  to  assess  as  to  magnitude  or  sign.  The  second  order 
terms  are  those  widely  used  in  approximate  theories  of  buckling,  but  are 
clearly  valid  only  in  a range  of  "moderately  large"  deflections.  At  final 
deflections  approaching  a third  of  the  span,  this  representation  is  quest- 
ionable. The  treatment  of  axial  forces  as  reactions  which  do  not  contri- 
bute to  the  energy  dissipation  rate  implies  a certain  over-estimate  of 
local  strength.  However  it  does  not  seem  possible  from  this  to  argue  how 
the  final  deflections  will  be  affected.  Intuitively  we  feel  that  the  er- 
rors due  to  these  aspects  are  not  serious  ones,  but  further  investigation 
is  necessary. 

(6)  The  pressure  pulse  is  idealized  as  "impulsive",  delivering  finite  impulse 
in  vanishingly  small  time.  It  is  well  known  that  this  idealization  leads 
to  over-estimates  of  major  deflections  [15]  when  applied  to  rigid-perfectly 
plastic  structures  using  small  deflection  equations.  The  error  in  the 
over-estimate  in  those  cases  can  be  shown  to  be  given  approximately 
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by  the  ratio  of  the  pulse  duration  to  the  duration  of  motion  of  the  struc- 
ture. These  results  certainly  do  not  apply  exactly  to  problems  where  visco 
plastic  and  large  deflection  effects  are  important.  However  in  the  present 
problems  the  pulse  duration  is  very  short  compared  to  the  response  time  of 
the  structure  (roughly  10  ysec  compared  to  10  msec),  and  the  resulting  er- 
ror is  believed  to  be  positive  but  of  minor  importance. 

(7)  The  integration  technique  used  to  obtain  final  deflections  of  the  structure 
with  initial  mode  form  velocity  employed  "instantaneous  mode"  solutions 

at  a sequence  of  times,  each  such  solution  corresponding  to  the  current  de- 
flection field.  They  are  linked  by  proper  differential  relations  only  at 
one  point  of  the  structure,  where  the  major  velocity  and  deflection  magni- 
tudes occur.  Hence  although  current  field  equations  are  satisfied,  some 
continuity  conditions  are  disregardsd.  Extremal  theorems  are  available  for 
the  small-small  equations  C7,  163.  These  can  be  applied  to  the  current 
state  and  show  that  the  mode  form  velocity  field  renders  the  energy  dis- 
sipation rate  an  extremum  among  all  kinematically  admissible  fields  with 
the  same  kinetic  energy.  Here  the  extremum  is  a minimum,  and  the  theorem 
characterizes  the  mode  form  solution  as  minimizing  the  rate  of  decrease  of 
kinetic  energy,  for  a fixed  level  of  kinetic  energy.  The  error  in  using 
this  technique  is  therefore  inferred  to  be  positive. 

(8)  Apart  from  the  basic  approximation  in  treating  the  actual  motion  as  a se- 
quence of  instantaneous  modes,  the  determination  of  current  mode  shapes 
and  accelerations  involves  solving  a nonlinear  eigenproblem  by  numerical 
means.  This  was  done  by  an  iterative  scheme,  as  outlined  in  the  Appendix. 
The  numerical  scheme  made  U3e  of  a subdivision  of  the  column  and  half-span 
into  100  elements.  It  furnishes  an  essentially  exact  solution  at  low  cost. 
No  difficulties  with  convergence  were  met  at  any  deflection  magnitudes,  and 
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it  is  believed  that  negligible  errors  are  introduced  in  this  part  of  the 
response  calculation. 

The  simplifications  discussed  above  are  involved  in  both  the  deflection 
bound  and  the  mode  technique,  with  the  exception  of  (7)  which  obviously  refers 
only  to  the  latter  approach.  It  may  be  noted  that  in  the  numerical  evaluation 
of  the  deflection  bound,  the  iterative  scheme  used  for  determining  the  deflected 
shape  under  static  loading  has  common  ground  with  that  used  in  the  numerical 
solution  of  the  eigen-problem  of  the  extended  mode  technique. 

It  is  seen  that  the  errors  that  can  be  identified  with  the  idealizations 
and  approximations  as  listed  are  mainly  positive,  i.e.,  such  as  to  cause  too 
large  a deflection.  The  exceptional  cases  are  items  (5)  and  (8),  where  no  state- 
ments on  the  sign  of  the  error  seem  tenable.  In  some  cases,  for  example  the 
neglect  of  elastic  deformation  and  of  strain  rate  history  effects , the  arguments 
are  somewhat  conjectural.  This  seems  unavoidable  at  the  present  stage. 

The  comparison  of  test  results  with  the  large  deflection  bounds  indicates 
that  the  method  of  calculation  of  the  bounds  is  satisfactory.  However,  their 
comparison  with  large  deflections  predicted  by  the  mode  technique  shows  incon- 
sistencies. In  the  tests  with  uniform  impulse  where  the  intrinsic  error  is 
positive,  as  are  most  of  the  other  errors,  the  estimated  deflection  curve  should 
lie  distinctly  above  the  test  points;  the  actual  test  points  lie  essentially 
on  the  curves.  Figs.  10,  11.  In  the  tests  with  concentrated  impulse,  the  in- 
trinsic error  is  negative.  This  is  shown  clearly  by  the  results  for  steel 
frames  in  Fig.  8,  which  suggests  that  the  other  errors  are  small  by  comparison. 
However,  the  results  for  the  titanium  frames  in  Fig.  9 are  then  anomalous, 
since  to  explain  them  the  incidental  errors  must  be  large.  Here  perhaps  the 
influence  of  elastic  deformations  is  more  important  as  indicated  by  the  energy 
ratio  criterion. 
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We  note  the  comparison  between  test  and  estimated  results  for  a secondary 
deflection  magnitude,  namely  the  inward  deflection  at  the  top  of  each  column. 
Figure  12  illustrates  these  results  for  the  titanium  frames  of  both  types  [1]. 

The  test  points  agree  reasonably  well  with  the  predictions  of  the  extended 
mode  technique.  In  a "small  deflection"  theory  this  deflection,  of  course, 
is  zero. 

Finally,  Fig.  13  illustrates  the  comparison  of  the  estimated  response  time 
t*  by  the  extended  mode  technique  with  the  two  times  obtained  from  strain  gage 
records  [1].  The  time  t£  is  the  time  for  the  first  maximum  strain.  The  time 
tj{  is  the  intercept  on  the  rising  strain  curve  of  a line  drawn  at  the  estimated 
final  (permanent)  strain.  The  estimated  response  time  falls  between  these  two 
measured  times  with  surprising  consistency. 

In  conclusion,  the  present  work  has  illustrated  the  application  to  i arge 
deflections  of  an  impulsively  loaded  viscoplastic  structure  of  bound  and  mode 
estimation  techniques.  The  efficiency  of  both  techniques  compared  to  purely 
numerical  approaches  has  been  emphasized.  Comparisons  with  fairly  comprehensive 
experimental  results  [l]  have  validated  the  deflection  bound  calculation,  but 
have  shown  up  inconsistencies  in  the  estimates  of  final  deflections  by  the  ex- 
tended mode  technique.  These  demonstrate  the  need  for  further  study  of  the 
approximations  made  in  implementing  this  method.  Despite  this,  it  is  worth 
emphasizing  that  the  estimated  deflections  by  the  extended  mode  method  in  most 
cases  agree  extremely  well  with  the  final  deflections  observed  in  the  tests. 

In  the  one  test  series  where  the  discrepancies  are  large  (Fig.  9),  the  estimates 
are  conservative  (actual  permanent  deflections  are  smaller);  these  occur  in  cir- 
cumstances of  relatively  small  final  deflections  where  elastic  deflections  are  more 
important . Hence  there  seems  no  doubt  that  practical  use  can  be  made  of  this 
approach,  with  normal  caution. 
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Captions  for  Fipures 


Frame  tvpes  considered  here  and  in  experiments  [1"].  Type  (a)  has 
steel  block  at  midpoint  of  transverse  member,  tvr>e  (h)  has  uniform 
transverse  member. 

Notation  for  mode  form  solution. 

Final  deflection  as  function  of  initial  mode  form  velocity  (both 
dimensionless):  (a)  showing  dependence  on  parameter  a,  (b)  show- 
ing dependence  on  power  n of  strain  rate  behavior  law. 

Time  of  response  as  function  of  initial  mode  form  velocity  (both 
dimensionless);  (a)  showing  dependence  on  parameter  a,  (b)  show- 
ing dependence  on  power  n of  strain  rate  behavior  law. 

Tvpical  curves  of  mode  velocity  as  function  of  time. 

"Master  response  curves"  of  final  deflection-thickness  ratio  as 
function  of  initial  mode  velocity  amplitude,  for  tvoe  (b)  frames. 

"Master  response  curves"  of  final  deflection-thickness  ratio  as 
function  of  initial  mode  velocity  amplitude,  *or  type  (a)  frames. 

Comparison  of  observed  final  midpoint  deflections  in  tests  [1] 
with  deflection  bound  and  predictions  of  mode  approximation  tech- 
nique, for  steel  frames  of  type  (a). 

Comparison  of  observed  final  midpoint  deflections  in  tests  fl] 
with  deflection  bound  and  predictions  of  mode  approximation  tech- 
nique, for  titanium  frames  of  type  (a). 

Comparison  of  observed  final  midpoint  deflections  in  tests  [1] 
with  deflection  bound  and  predictions  of  mode  approximation  tech- 
nique, for  steel  frames  of  type  (b). 

Comparison  of  observed  final  midpoint  deflections  in  tests  [1] 
with  deflection  bound  and  predictions  of  mode  approximation  tech- 
nique, for  titanium  frames  of  type  (b). 

Comparison  of  final  lateral  deflection  at  top  of  column  observed 
in  tests  [1]  with  predictions  of  the  extended  mode  approximation 
technique . 

Comparison  of  time  Lo  final  deflection  predicted  by  extended  mode 
techninue  with  times  obtained  from  strain  jjape  records  in  experi- 
ments Cl];  t^  is  time  of  peak  deflection  (strain),  t^  is  time 
of  intercept  with  rising  strain-time  curve  of  straight  line  drawn 
at  final  strain  level;  strain  measured  at  ton  of  column. 
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Appendix 

The  iterative  scheme  used  to  determine  the  instantaneous  mode  at  each 
time  stage  in  the  mode  technique  is  essentially  the  same  as  that  used  in 
the  somtion  for  finite  deflections  of  a circular  plate  T4].  It  is  note- 
worthy that  the  deflection  function  giving  the  deflection  bound  can  be 
obtained  by  the  same  iteration  scheme.  We  shall  sketch  briefly  how  the  equa- 
tions of  the  two  problems  can  be  put  in  analogous  form,  omitting  most  of 
the  details. 

Equations  only  for  the  transverse  member  BC  will  be  written,  for 
brevity.  We  look  first  at  the  mode  technique.  When  the  modal  forms  of 
Eqs.  (20)  are  used,  the  dimensionless  bending  moment  in  the  transverse 
member  BC  can  be  written  as 

1/n' 

m = w(-‘  ) | -4>”  | spn(  -<j>”)  (Al) 

Equations  (19b),  (15d),  (15b)  can  be  integrated  with  respect  to  x and 
put  in  the  following  forms: 


(19b) 

(15d) 

(15b) 


s = <02dx  + sq] 

1 Jo 

1/n'  1/n'  , 

p(—  ) [ sgn(-if)^)]  + 4sw ' 


(A2) 

(A3) 


••  fX 

w.J  <t>  Jx  + A] 
o 

(A4 ) 


Here  sq  and  A are  constants  and  the  prime  denotes  differentiation  with 
respect  to  x,  position  in  member  BC  measured  from  B.  (All  quantities  refer 
to  this  member,  but  ripts  are  omitted).  Similar  equations  can  be 

written  for  member  AB,  and  the  boundary  conditions  include  fixing  at  A and 


...  B c | 
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C,  continuity  at  B , and  the  normalizing  condition  ^ = 1 at  C . After 
use  of  Eq.  (A3),  Eq.  (A4)  can  be  written  as 


1/n'  , rrx  r (X 

sgn(-4>” )( | -«|>^ | ) = Aj  (J^dx  + A - 4w'  Ur-  <J> 

[_  1 

1/n' 

# • •*/»(?) 


where 


2dx  * So 


Integrating  Eq.  (A2),  we  obtain 


*2  = - ^ [|oW’ndx  - |oW'*i  dx]  * 


At  any  stage  of  the  response  the  slope  function  w'(x)  can  be 
estimated  from  the  field  at  a preceding  stage  together  with  the  velocity  and 
acceleration  fields  at  that  stage.  Admissible  shape  functions  ^(x)  in 
BC  and  $2(y)  in  AB  can  also  be  guessed,  satisfying  the  stated  conditions 
at  A,  B and  C;  the  functions  for  the  preceding  stage  may  be  used,  except 
at  t = 0 . Eq.  (A6)  then  furnishes  <J>2(x)  in  BC  » apart  from  integration 
constants.  With  use  of  the  analogous  equations  for  AB  and  the  fixing, 
continuity,  and  normalization  conditions,  all  of  the  integration  constants 
can  be  determined,  and  new  shape  functions  obtained.  From  each  set  of 
functions  the  quantity  A (Eq.  ( A5 ) ) can  be  computed  from  the  energy 
rate-dissipation  equation  (14),  which  takes  the  form 


f1  Lo  r1 

j |-^|  n*  dx  + t—  f |-*^|  n'  dy 
0 1 ' o 

J1  2 2 L2  f1  2 2 

(^  + 4>2 )dx  + — J («(>1  + <f>2)dy  + 
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Cycles  of  iteration  may  be  terminated  when  the  magnitude  A reaches  a 
steady  value.  From  the  final  values,  the  velocity  and  displacement  fields 


can  be  written,  and  approximate  functions  for  a subsequent  instant 
obtained. 

At  the  start,  w'  = 0 and  is  chosen  arbitrarily,  with  guessed 

admissible  functions  and  $^(x)  in  AB  and  BC  , respectively. 

The  response  terminates  at  ^(t^)  = 0 , final  deflections  being  obtained 

by  interpolation  when  wA  first  becomes  negative. 

The  deflection  bound  requires  only  one  determination  of  shape  functions 

over  the  frame,  rather  than  functions  at  a sequence  of  times  during  the 

response  as  in  the  mode  technique.  To  obtain  an  upper  bound  on  the  final 

deflection  of  the  midpoint  C of  the  frame,  a force  P is  applied  at 

this  point.  If  the  work  done  quasi-statically  by  this  force  in  time  tf 

is  not  less  than  K , where  t^  is  the  actual  response  rime  and  K the  given 
o f o 

initial  kinetic  energy,  then  the  deflection  corresponding  to  this  force  is 
an  upper  bound  on  the  same  component  of  deflection  resulting  from  the  given 
initial  velocities. 

Use  of  minimum  work  paths  to  fixed  final  strains  at  time  tf  [2,3] 
for  the  material  behavior  represented  by  Eqs.  (11,12)  enables  us  to  write 
moment -curvature  relations 


Kt-  = < = x t _ 

f or 


n'+l  uM 


sgn  M 


(A8) 


where  < f and  Mr  are  curvature  and  moment  at  time  t^  , respectively. 
Nondimensional  stresses  may  be  defined  as 
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where  t = 2L, /p/o 
1 0 


is  the  reference  time  used  also  in  the  mode  method. 
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The  equation  of  equilibrium  for  the  beam  member  BC  with  a concentrated 
force  P at  L^,  written  as  a delta  function,  is 

M"  ♦ (Nw>'  = - . (A10) 

This  may  be  integrated  once,  and  with  nondimens ional  variables  as  in  Eqs. 
(9)  and  (A9)  the  four  equations  for  this  member  can  be  written  as 


m ' + 4sw ' 
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( Alld ) 


where  p = 


1/n'  L,P 


-i-  , and  a = 8L,e  /p/o  /H2  . The  analogous 

2M  * loo 


equations  for  AC  , and  the  boundary  and  continuity  equations,  are  omitted 
for  brevity.  The  work-kinetic  energy  condition  Eq.  (lb)  takes  the  form 


ST  5 2w°  ^ 


»l/n 


[ I1  |i|"’wdx  t £ | 1 l|r,Hhldy 

•O  1 • 0 


(A12) 


In  these  equations,  the  response  time  t^  appears  explicitly  only  in 
Eq.  (A12),  because  of  the  definitions  of  m and  p . The  inequality  shows 
that  tf  , which  is  unknown,  may  be  replaced  by  an  upper  bound  tjj  . We 
shall  not  attempt  to  derive  an  exact  upper  bound.  Note  that  t^  appears 
only  as  t^/n'  , where  n'  is  large  (e.g.,  10),  so  results  are  insensitive 
to  tf  . For  small  deflections  and  rigid-perfectly  plastic  behavior,  upper 
bounds  on  the  response  time  for  the  present  frame  problems  can  be  obtained 


by  Lee's  method  [17]  and  written  as 
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where  n = 2k/l+l/3k  for  type  (a)  frames  (attached  mass  and  concentrated 
impulse)  and  n = 2//T  for  type  (b)  frames.  The  effect  of  finite  deflections 
is  to  lengthen,  while  that  of  strain  rate  sensitivity  is  to  shorten  the 
response  time.  Hence  while  the  expressions  for  t^  are  not  exact  for 
our  conditions,  they  are  probably  more  than  adequate. 

It  is  convenient  to  write  the  deflection  field  in  the  form 

w = D<fr^(x)  in  BC  ; w = D^(y)  in  AB  (A14a) 

u = D$2(x)  in  BC  ; u = D$2(y)  in  AB  (A14b) 

where  D is  the  dimensionless  displacement  at  the  midpoint  C of  the 
beam  member.  An  iteration  scheme  is  started  by  giving  D an  arbitrary 
value,  e.g.,  D = 20  , and  assuming  numerical  functions  $^(x)  in  BC  and 
♦2(y)  in  AB  . Equations  such  as  Eq.  (Alld)  then  furnish  the  correspond- 
ing axial  components  $2(x)  in  BC  and  $^(y)  in  AB  . The  equations 
of  equilibrium  (e.g.,  Eq.  (Alla)),  together  with  the  moment-curvature 
equations  (e.g.,  Eq.  (Allb)),  can  now  be  written  entirely  in  terms  of  the 

transverse  shape  functions  ^(x)  and  2 ( y ) and  their  derivatives, 

so  that  from  starting  functions  new  functions  are  obtained  by  quadratures, 
with  constants  of  integration  which  must  be  found  from  the  boundary  and 
continuity  equations.  The  forms  of  the  equations  are  identical  with  those 
in  the  mode  method,  with  different  meaning  of  the  constants. 

The  experiments  with  which  we  have  made  comparisons  involve  direct 
measurements  of  the  total  impulse  (I)  applied  to  the  frame.  Thus  for  the 
two  frame  types 


Using  these  and  the  expressions  for  the  upper  bound  on  response  time,  Eq.  (A13), 
we  have  relations  between  impulse  and  midpoint  deflection  ratio  D = w^/ri 
as  follows 
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where  (for  both  types) 
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Suitable  values  of  y = o'/o  and  v = n'/n  are  obtained  from  the  initial 

o o 

strain  rates  by  means  of  Eq.(ll).  The  initial  mode  form  velocity  field,  with 
•o 

magnitude  wft  derived  from  given  initial  velocity  or  impulse  values  through 
the  mode  matching  technique  of  Eq.(6),  may  be  used  with  $(x)  and  +(y)  determined 
as  outlined  above.  This  field  gives  initial  maximum  curvature  rates  and  strain 
rates  which,  of  course,  are  not  the  actual  initial  strain  rates,  but  are 
smoothed  values  that  are  appropriate  for  present  purposes. 
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